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Abstract 

We present a preliminary measurement of time-dependent CP-violating asymmetries in B° — > 
J/t/jK® and B° — ► ip(2S)Kg decays recorded by the BABAR detector at the PEP-II asymmetric B 
Factory at SLAC. The data sample consists of 9.0 fb _1 collected at the T(4S) resonance and 0.8 fb _1 
off-resonance. One of the neutral B mesons, produced in pairs at the T(4S), is fully reconstructed. 
The flavor of the other neutral B meson is tagged at the time of its decay, mainly with the charge 
of identified leptons and kaons. The time difference between the decays is determined by measuring 
the distance between the decay vertices. Wrong-tag probabilities and the time resolution function 
are measured with samples of fully-reconstructed semileptonic and hadronic neutral B final states. 
The value of the asymmetry amplitude, sin2/3, is determined from a maximum likelihood fit to the 
time distribution of 120 tagged B° J/ipK® and B° -> ^(2S)K Q S candidates: sin2/3=0.12 ± 0.37 
(stat) ± 0.09 (syst). 
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1 Introduction 

The BABAR detector at the PEP-II asymmet- 
ric B Factory at SLAC has been taking data 
since the end of May, 1999. As of the date of 
this conference a data sample of ~9.0fb~ 1 has 
been collected at the T(4S) resonance, with an 
additional 0.8 fb -1 taken off-resonance. A vari- 
ety of preliminary B physics results using this 
data sample have been submitted to this con- 
ference and reported in parallel sessions. This 
paper will discuss in some detail BABAR's first 
measurements of CP-violating asymmetries in 
B° -► J/4)K° S and B° -> 4>{2S)K° decays. 

2 Motivation and Overview 

The CP- violating phase of the three-genera- 
tion Cabbibo-Kobayashi-Maskawa 
(CKM) quark mixing matrix can provide an el- 
egant explanation of the well-established CP- 
violating effects seen in decay However, 
studies of CP violation in neutral kaon decays 
and the resulting experimental constraints on 
the parameters of the CKM matrix j|] do not, 
in fact, yet provide a test of whether the CKM 
phase describes CP violation Q. 

The unitarity of the three generation CKM 
matrix can be expressed in geometric form as six 
triangles of equal area in the complex plane. A 
nonzero area[Q] directly implies the existence of 
a CP-violating CKM phase. The most experi- 
mentally accessible of the unitarity relations, in- 
volving the two smallest elements of the CKM 
matrix, V u b and Vtd, has come to be known as 
the (B) Unitarity Triangle. Because the lengths 
of the sides of this Unitarity Triangle are of the 
same order, the angles can be large, leading to 
potentially large CP-violating asymmetries from 
phases between CKM matrix elements. 

The CP- violating asymmetry in b — > tics de- 
cays of the B° meson such as B°/B° — ► J/ip K® 
(or B°/B° -> vp(2S)K°) is caused by the inter- 
ference between mixed and unmixed decay am- 
plitudes. A state initially prepared as a B° (B ) 
can decay directly to J/tp K® or can oscillate into 



a B° (B°) and then decay to J/ipK° s . With 
little theoretical uncertainty, the phase differ- 
ence between these amplitudes is equal to twice 
the angle (5 = &ig[-V cd V; b /V td V t t] of the Uni- 
tarity Triangle. The CP-violating asymmetry 
can thus provide a crucial test of the Standard 
Model. The interference between the two ampli- 
tudes, and hence the CP asymmetry, is maximal 
when the mixing probability is at its highest, 
i.e., when the lifetime t is approximately 2.2 B° 
proper lifetimes. 

In e + e~ storage rings operating at the T(4S) 
resonance a B°B° pair produced in T(4S) de- 
cay evolves in a coherent P-wave until one of 
the B mesons decays. If one of the B mesons 
(Btag) can be ascertained to decay to a state of 
known flavor at a certain time U a g, the other B 
is at that time known to be of the opposite fla- 
vor. For this measurement, the other B (Bcp) 
is fully reconstructed in a CP eigenstate (J/ipK° 
or ip(2S)Kg). By measuring the proper time in- 
terval At = tcp — ttag from the B tag decay time 
to the decay of the Bcp, it is possible to deter- 
mine the time evolution of the initially pure B° 
or B° state. The time-dependent rate of decay 
of the Bcp final state is given by 

/±(Ai; T, Am d , V sin 2(5) = ir e - r l A 'l 

x [ 1 ± V sin 2(5 x sin Am d At\ , (1) 

where the + or — sign indicates whether the 
B tag is tagged as a B° or a B°, respectively. 
The dilution factor T> is given by T> = 1 — 2w, 
where w is the mistag fraction, i.e., the probabil- 
ity that the flavor of the tagging B is identified 
incorrectly. A term proportional to cos Am d At 
would arise from the interference between two 
decay mechanisms with different weak phases. 
In the Standard Model, the dominant diagrams 
(tree and penguin) for the decay modes we con- 
sider have no relative weak phase, so no such 
term is expected. 

To account for the finite resolution of the de- 
tector, the time-dependent distributions f± for 
B° and B° tagged events (Eq. |l|) must be convo- 
luted with a time resolution function lZ(At; a): 
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F±( At ; T, Am d , V sin 2/3, a) = 

f±(At;T, Am d ,V sin 2/3) (8) TZ( At; a) , (2) 

where a represents the set of parameters that 
describe the resolution function. 

In practice, events are separated into differ- 
ent tagging categories, each of which has a dif- 
ferent mean dilution T>i , determined individually 
for each category. 

It is possible to construct a CP- violating ob- 
servable 



Acp(At) 



F+{At) - T-(At) 
F+(At) + F-{At) 



which is proportional to sin2/3: 

Acp{At) ~ V sin 2\3 x sinAm d At 



(3) 



(4) 



Since no time-integrated CP asymmetry 
effect is expected, an analysis of the time- 
dependent asymmetry is necessary. At an 
asymmetric-energy B Factory, the proper decay- 
time difference At is, to an excellent approxima- 
tion, proportional to the distance Az between 
the two _B°-decay vertices along the axis of the 
boost, At Az/c{Pi). At PEP-II the average 
boost of B mesons, (/?7), is 0.56. The distance 
Az is 250 ^m per B° lifetime, while the typical 
Az resolution for the BABAR detector is about 
110 fim. 

Since the amplitude of the time-dependent 
CP-violating asymmetry in Eq. |] involves the 
product of V and sin2/3, one needs to deter- 
mine the dilution factors T>i (or equivalently 
the mistag fractions tOj) in order to extract the 
value of sin2/3. The mistag fractions can be 
extracted from the data by studying the time- 
dependent rate of B°B° oscillations in events in 
which one of the neutral B mesons is fully recon- 
structed in a self-tagging mode and the other B 
(the B tag ) is flavor-tagged using the standard 
CP analysis flavor-tagging algorithm. In the 
limit of perfect determination of the flavor of 
the fully-reconstructed neutral B, the dilution in 
the mixed and unmixed amplitudes arises solely 



from the Btag side, allowing the values of the 
mistag fractions Wi to be determined. 

The value of sin2/3 is extracted by maximiz- 
ing the likelihood function The value of sin2/3 is 
extracted by maximizing the likelihood function 

ln£cp = 



Ej E BQtag ln F+ ( At 5 r > Am d, a, T>i sin 2/3) 
+ Eso t ln ^"-(At; T, Am d , a, T>i sin 2(3 ) 



where the outer summation is over tagging cat- 
egories i. 

2.1 Overview of the analysis 

The measurement of the CP-violating asymme- 
try has five main components : 

• Selection of the signal B°/B° -> J/ipK Q s 
and B°/B° -» ip(2S)K° events, as de- 
scribed in detail in Ref. 5. 

• Measurement of the distance Az between 
the two B° decay vertices along the T(4S) 
boost axis, as described in detail in Refs. 6 
and 7. 



Determination of the flavor of the Bt ag , as 
described in detail in Ref. 6. 



• Measurement of the dilution factors T>i 
from the data for the different tagging cat- 
egories, as described in detail in Ref. 6. 

• Extraction of the amplitude of the CP 
asymmetry and the value of sin2/3 with an 
unbinned maximum likelihood fit. 

Whenever possible, we determine time and 
mass resolutions, efficiencies and mistag frac- 
tions from the data. 

3 Sample selection 

For this analysis we use a sample of 9.8 fb _1 of 
data recorded by the BABAR detector^] between 
January 2000 and the beginning of July 2000, of 



(5) 
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which 0.8 fb was recorded 40MeV below the 
T(4:S) resonance (off-resonance data). 

A brief description of the BABAR detector 
and the definition of many general analysis pro- 
cedures can be found in Ref. 8. Charged parti- 
cles are detected and their momenta measured 
by a combination of a central drift chamber 
(DCH) filled with a helium-based gas and a 
five-layer, doubled-sided silicon vertex tracker 
(SVT), in a 1.5 T solenoidal field produced 
by a superconducting magnet. The charged 
particle momentum resolution is approximately 
{5pt/pt) 2 = (0.0015 p T ) 2 + (0.005) 2 , where p T 
is measured in GeV/c. The SVT, with typical 
10 /urn single-hit resolution, provides vertex in- 
formation in both the transverse plane and in 
the z direction. Vertex resolution is typically 
50 fj,va in z for a fully reconstructed B meson, 
depending on the decay mode, and of order 100 
to 150 /mi for a generic B decay. Leptons and 
hadrons are identified with measurements from 
all the BABAR components, including the energy 
loss dE/dx from a truncated mean of up to 40 
samples in the DCH and at least 8 samples in 
the SVT. Electrons and photons are identified 
in the barrel and the forward regions by the Csl 
electromagnetic calorimeter (EMC). Muons are 
identified in the instrumented flux return (IFR) . 
In the central polar region the Cherenkov ring 
imaging detector (DIRC) provides K-ir separa- 
tion with a significance of at least three standard 
deviations over the full momentum range for B 
decay products above 250 MeV/ c. 

3.1 Particle identification 

An electron candidate must be matched to an 
electromagnetic cluster of at least three crys- 
tals in the Csl calorimeter. The ratio of the 
cluster energy to the track momentum, E/p, 
must be between 0.88 and 1.3. The lateral mo- 
ment of the cluster must be between 0.1 and 0.6, 
and the Zernike moment of order (4,2)Qmust be 
smaller than 0.1. In addition the electron can- 

lr The lateral and Zernike moments are cluster shape 
variables introduced in Ref. 8. 



didate track in the drift chamber must have a 
dE/dx measurement consistent with that of an 
electron and, if measured, the Cherenkov angle 
in the DIRC must be consistent with that of an 
ultra-relativistic particle. 

Muon identification relies principally on the 
measured number of interaction lengths, N\, 
penetrated by the candidate in the IFR iron, 
which must have a minimum value of 2.2 and, at 
higher momenta, must be larger than jVf — 1, 
where N^ xp is the expected number of interac- 
tion lengths for a muon. The number of IFR 
layers with a "hit" must be larger than two. To 
reject hadronic showers, we impose criteria on 
the number of IFR strips with a hit as a function 
of the penetration length, and on the distance 
between the strips with hits and the extrapo- 
lated track. In the forward region, which suf- 
fers from accelerator-related background, extra 
hit-continuity criteria are applied. In addition, 
if the muon candidate is in the angular region 
covered by the EMC, the energy deposited by 
the candidate in the calorimeter must be larger 
than 50 MeV and smaller than 400 MeV. (The 
expected energy deposited by a minimum ioniz- 
ing particle is about 180 MeV.) 

Particles are identified as kaons if the ratio 
of the combined kaon likelihood to the combined 
pion likelihood is greater than 15. The com- 
bined likelihoods are the product of the indi- 
vidual likelihoods in the SVT, DCH and DIRC 
subsystems. In the SVT and DCH tracking de- 
tectors, the likelihoods are based on the mea- 
sured dE/dx truncated mean compared to the 
expected mean for the K and tt hypotheses, 
with an assumed Gaussian distribution. The 
dE/dx resolution is estimated on a track- by- 
track basis, based on the direction and momen- 
tum of the track and the number of energy de- 
position samples. For the DIRC, the likelihood 
is computed by combining the likelihood of the 
measured Cherenkov angle compared to the ex- 
pected Cherenkov angle for a given hypothesis, 
with the Poisson probability of the number of 
observed Cherenkov photons, given the number 
of expected photons for the same hypothesis. 
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DIRC information is not required for particles 
with momentum less than 0.7GeV/c, where the 
DCH dE/dx alone provides good K/ir discrim- 
ination. 

3.2 Data samples 

We define three event classes^: 

• the CP sample, containing B° candi- 
dates reconstructed in the CP eigenstates 
J/ipK® or ip(2S)K®. The charmonium 
mesons J/0 and ip(2S) are reconstructed 
through their decays to e + e~ and /i + pT . 
The tp(2S) is also reconstructed through 
its decay to J/07r + 7r~. The K® is recon- 
structed through its decay to 7r + ir~ and 
7r°vr . The selection criteria for the CP 
sample are described in the next section. 

• the fully reconstructed B° samples, con- 
taining B° candidates in either semilep- 
tonic or hadronic flavor eigenstates. The 
sample of semileptonic decays contains 
candidates in the B — > D*~i + vi mode 
{i + = e + or the sample of hadronic 
neutral decays contains B° candidates 
in the D^~tt + , D^~p+ and D^~af 
modes; the sample of hadronic charged de- 
cays contains B + candidates in the D n + , 
and Z)*°7r+ (with D*° -► ir°D°) modes. 
The selection criteria for these samples 
are described in Refs. 6 and 7.We recon- 
struct 7500 B° D*~t + vi candidates, 
~ 2500 candidates in hadronic B° final 
states, and ss 2300 candidates in hadronic 
B + final states. 

• the charmonium control samples, contain- 
ing fully reconstructed neutral or charged 
B candidates in two-body decay modes 
with a J/0 in the final state, such as 
B + -> J/ipK+ or B° -> J/ip(K*° -» 
K + ir~). The selection criteria for these 

2 Throughout this paper, conjugates of fiavor- 
eigenstate modes are implied. 



samples are described in Ref. 5. We recon- 
struct 570 B + — > J/0 K + candidates and 
237 B° -» J/4j(K*° -» K+7T-) candidates. 

Signal event yields and purities for the indi- 
vidual samples are summarized in Table |l]. 

3.3 The CP sample 

We select events with a minimum of four re- 
constructed charged tracks in the region defined 
by 0.41 < 6 i a b < 2.41. Events are required 
to have a reconstructed vertex within 0.5 cm 
of the average position of the interaction point 
in the plane transverse to the beamline, and 
a total energy greater than 5 GeV in the fidu- 
cial regions for charged tracks and neutral clus- 
ters. To reduce continuum background, we re- 
quire the second-order normalized Fox- Wolfram 
moment H (i?2 = #2/ '-Ho) of the event to be less 
than 0.5. 

The selection criteria for the J/ipK® and 
ip(2S)K® events are optimized by maximizing 
the ratio S /y/S + £>, where S (the number of 
signal events that pass the selection) is deter- 
mined from signal Monte Carlo events, and B 
(the number of background events that pass 
the selection) is estimated from a luminosity- 
weighted average of continuum data events and 
nonsignal BB Monte Carlo events. 

For the J/0 or tp(2S) — > e + e~ candidates, 
at least one of the decay products is required 
to be positively identified as an electron or, if 
outside the acceptance of the calorimeter, to 
be consistent with an electron according to the 
drift chamber dE/dx information. If both tracks 
are within the calorimeter acceptance and have 
a value of E/p larger than 0.5, an algorithm 
for the recovery of Bremsstrahlung photons || 
is used. 

For the J/0 or tp(2S) — ► candidates, 
at least one of the decay products is required to 
be positively identified as a muon and the other, 
if within the acceptance of the calorimeter, is re- 
quired to be consistent with a minimum ionizing 
particle. 

We select J/0 candidates with an invariant 
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mass greater than 2.95GeV/c 2 and 3.06GeV/c 2 
for the e + e~ and /J + /U~ modes, respectively, and 
smaller than 3.14GeV/c 2 in both cases. The 
ip(2S) candidates in leptonic modes must have a 
mass within 50MeV/c 2 of the ijj(2S) mass. The 
lower bound is relaxed to 250MeV/c 2 for the 
e + e~ mode. 

For the vp(2S) — > J/ipir + ir~ mode, mass- 
constrained J/0 candidates are combined with 
pairs of oppositely charged tracks considered as 
pions, and tp(2S) candidates with mass between 
3.0GeV/c 2 and 4.1GeV/c 2 are retained. The 
mass difference between the ip(2S) candidate 
and the J/0 candidate is required to be within 
15MeV/c 2 of the known mass difference. 

Kg candidates reconstructed in the 7r + 7r~ 
mode are required to have an invariant mass, 
computed at the vertex of the two tracks, 
between 486MeV/c 2 and 510MeV/c 2 for the 
J/0 Kg selection, and between 491MeV/c 2 and 
505MeV/c 2 for the ^{2S)K° S selection. 

For the JfyK® mode, we also consider the 
decay of the Kg into 7r°7r°. Pairs of ir° candi- 
dates, with total energy above 800 MeV and in- 
variant mass, measured at the primary vertex, 
between 300 and 700 MeV/ c 2 , are considered as 
Kg candidates. For each candidate, we deter- 
mine the most probable Kg decay point along 
the path defined by the Kg momentum vector 
and the primary vertex of the event. The decay- 
point probability is the product of the \ 2 prob- 
abilities for each photon pair constrained to the 
7T° mass. We require the distance from the de- 
cay point to the primary vertex to be between 
— 10 cm and +40 cm and the K° mass measured 
at this point to be between 470 and 536 MeV/ c 2 . 

Bqp candidates are formed by combining 
mass-constrained J/0 or ip(2S) candidates with 
mass-constrained Kg candidates. The cosine of 
the angle between the Kg three-momentum vec- 
tor and the vector that links the J/0 and Kg 
vertices must be positive. The cosine of the he- 
licity angle of the J/0 in the B rest frame must 
be less than 0.8 for the e + e~ mode and 0.9 for 
the fi + fi~ mode. 

For the t/j(2S)Kg candidates, the helicity 
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Beam-energy substituted mass (GeV/c~) 
Figure 1: J/tpK® (Kg — * 7r + 7r~) signal. 
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angle of the ip{2S) must be smaller than 0.9 
for both leptonic modes. The Kg flight length 
with respect to the i/>(2S) vertex is required 
to be greater than 1mm. In the ip(2S) — > 
J/ip Tr + ir~ mode, the absolute value of the cosine 
of the angle between the Bqp candidate three- 
momentum vector and the thrust vector of the 
rest of the event, in the center-of-mass frame, 
must be less than 0.9. 

Bcp candidates are identified with a pair 
of nearly uncorrelated kinematic variables: the 
difference AE between the energy of the Bqp 
candidate and the beam energy in the center-of- 
mass frame, and the beam-energy substituted 
massQ msE- The signal region is defined by 
5.270 GeV/c 2 < m SE < 5.290 GeV/c 2 and an 
approximately three-standard-deviation cut on 



AE (typically \AE\ < 35MeV). 

Distributions of AE and mgE are shown in 
Fig. Q Q and || for the CP samples and in Fig. |4| 
and H for the charmonium control samples. Sig- 
nal event yields and purities, determined from 
a fit to the msE distributions after selection on 
AE, are summarized in Table |l[ 

The CP sample used in this analysis is com- 
posed of 168 candidates: 121 in the J/ip 
(K° -» vr+vr-) channel, 19 in the J/^K Q S {K Q S -> 
vr vr ) channel and 28 in the ip(2S)K° s -> 
7r + 7r _ ) channel. 

4 Time resolution function 

The resolution of the At measurement is dom- 
inated by the z resolution of the tagging ver- 
tex. The tagging vertex is determined as fol- 
lows. The three- momentum of the tagging B 
and its associated error matrix are derived from 
the fully reconstructed Bqp candidate three mo- 
mentum, decay vertex and error matrix, and 
from the knowledge of the average position 
of the interaction point and the T(4S) four- 
momentum. This derived Bt ag three-momentum 
is fit to a common vertex with the remaining 
tracks in the event (excluding those from Bcp)- 
In order to reduce the bias due to long-lived 
particles, all reconstructed V° candidates are 
used as input to the fit in place of their daugh- 
ters. Any track whose contribution to the % 2 
is greater than 6 is removed from the fit. This 
procedure is iterated until there are no tracks 
contributing more than 6 to the x 2 or until 
all tracks are removed. Events are rejected if 
the fit does not converge for either the Bqp or 
B tag vertex. We also reject events with large 
Az ( \Az\ > 3mm) or a large error on Az 
(<7Az > 400 fxm). 

The time resolution function is described ac- 
curately by the sum of two Gaussian distribu- 
tions, which has five independent parameters: 



U( At; a) 

2 

E 



fi 



\ ai\f2ir 



exp(-(At-<5i) 2 /2ov 2 



(6) 
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Table 1: Event yields for the different samples used in this analysis, from the fit to ttt-se distribu- 
tions after selection on AE. The purity is quoted for ttise > 5.270 MeV/c 2 (except for D*~l + u). 



r\ l 

Sample 


t mal state 


Yield 


Purity (%) 


CP 


J/ipK% (Kg -► vr+Tr") 


124±12 


96 




J/fy-f^ -► 7r°7r°) 


18±4 


91 




v( 2 S) K s 


27±b 


93 


Hadronic 


D*~7T + 


622±27 


90 


(neutral) 


D*-p + 


419±25 


84 




J ~\ ;■■ -\- 

D a{ 


239±l9 


79 




U 71 


DoUznZD 


on 




u p 




84 




D~af 


225±20 


74 




total 


2438±57 


85 


Hadronic 




1755±47 


88 


(charged) 


D*TT+ 


543±27 


89 




total 


2293±54 


88 


Semileptonic 


D*-£ + u 


7517±104 


84 


Control 


j/i)K + 


597±25 


98 




^(2S)K + 


92±10 


93 




J/t/jK*° (K*° -► K+-K-) 


251±16 


95 



A fit to the time resolution function in Monte 
Carlo simulated events indicates that most of 
the events (/i = 1 — /2 = 70%) are in the core 
Gaussian, which has a width o\ « 0.6 ps. The 
wide Gaussian has a width 02 ~ 1.8 ps. Tracks 
from forward-going charm decays included in 
the reconstruction of the B tag vertex introduce a 
small bias, 8\ « — 0.2ps, for the core Gaussian. 

A small fraction of events have very large 
values of Az, mostly due to vertex reconstruc- 
tion problems. This is accounted for in the 
parametrization of the time resolution function 
by a very wide unbiased Gaussian with fixed 
width of 8ps. The fraction of events populat- 
ing this component of the resolution function, 
f w , is estimated from Monte Carlo simulation 
as ~ 1%. 

In likelihood fits, we use the error a At on 
At that is calculated from the fits to the two 
B vertices for each individual event. However, 
we introduce two scale factors Si and S2 for the 
width of the narrow and the wide Gaussian dis- 



tributions (oi = Si x UAt and 02 = Si x a^t) 
to account for the fact that the uncertainty on 
At is underestimated due to effects such as the 
inclusion of particles from D decays and possi- 
ble underestimation of the amount of material 
traversed by the particles. The scale factor Si 
and the bias 81 of the narrow Gaussian are free 
parameters in the fit. The scale factor £2 and 
the fraction of events in the wide Gaussian, /2, 
are fixed to the values estimated from Monte 
Carlo simulation by a fit to the pull distribution 
(<S 2 = 2.1 and f 2 = 0.25). The bias of the wide 
Gaussian, 62, is fixed at Ops. The remaining set 
of three parameters: 

a = {Si, Si, f w } (7) 

are determined from the observed vertex distri- 
bution in data. 

Because the time resolution is dominated by 
the precision of the B tag vertex position, we 
find no significant differences in the Monte Carlo 
simulation of the resolution function parameters 
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for the various fully reconstructed decay modes, 
validating our approach of determining the reso- 
lution function parameters a with the relatively 
high-statistics fully-reconstructed B° data sam- 
ples, and fixing these parameters in the likeli- 
hood fit for the determination of sin2/3 with the 
low-statistics CP sample. The differences in the 
resolution function parameters in the different 
tagging categories are also small. 

Table || presents the values of the parameters 
obtained from a fit to the hadronic B° sample. 
These values are used in the final fit for sin2/3. 

Table 2: Parameters of the resolution function 
determined from the sample of events with fully- 
reconstructed hadronic B candidates. 



Parameter 


Value 


Si (ps) 


-0.20 ±0.06 


from fit 




1.33 ±0.14 


from fit 


U (%) 


1.6 ±0.6 


from fit 


h (%) 


75 


fixed 


S 2 ( ps) 





fixed 


s 2 


2.1 


fixed 



5 B flavor tagging 

Each event with a CP candidate is assigned 
a B° or B° tag if the rest of the event (i.e., 
with the daughter tracks of the Bqp removed) 
satisfies the criteria for one of several tagging 
categories. The figure of merit for each tag- 
ging category is the effective tagging efficiency 
Qi = Si (1 — 2wi) 2 , where £j is the fraction of 
events assigned to category i and Wi is the prob- 
ability of misclassifying the tag B° or B° for 
this category. Wi is called the mistag fraction. 
The statistical error on sin2/3 is proportional to 
1/y/Q, where Q = J2iQi- 

Three tagging categories rely on the pres- 
ence of a fast lepton and/or one or more charged 
kaons in the event. Two categories, called neural 
network categories, are based upon the output 
value of a neural network algorithm applied to 



events that have not already been assigned to 
lepton or kaon tagging categories. 

In the following, the tag refers to the B tag 
candidate. In other words, a B° tag indicates 
that the Bcp candidate was in a B° state at 
At = 0; a B° tag indicates that the Bcp candi- 
date was in a B° state. 

5.1 Lepton and kaon tagging cate- 
gories 

The three lepton and kaon categories are called 
Electron, Muon and Kaon. This tagging tech- 
nique relies on the correlation between the 
charge of a primary lepton from a semileptonic 
decay or the charge of a kaon, and the flavor 
of the decaying b quark. A requirement on the 
center-of-mass momentum of the lepton reduces 
contamination from low-momentum opposite- 
sign leptons coming from charm semileptonic 
decays. No similar kinematic quantities can be 
used to discriminate against contamination from 
opposite-sign kaons. Therefore, for kaons the 
optimization of Q relies principally on the bal- 
ance between kaon identification efficiency and 
the purity of the kaon sample. 

The first two categories, Electron and Muon, 
require the presence of at least one identified 
lepton (electron or muon) with a center-of-mass 
momentum greater than l.lGeV/c. The mo- 
mentum cut rejects the bulk of wrong-sign lep- 
tons from charm semileptonic decays. The value 
is chosen to maximize the effective tagging effi- 
ciency Q. The tag is B° for a positively-charged 
lepton, B° for a negatively-charged lepton. 

If the event is not assigned to either the 
Electron or the Muon tagging categories, the 
event is assigned to the Kaon tagging category if 
the sum of the charges of all identified kaons in 
the event, T,Qk, is different from zero. The tag 
is B° if T,Qk is positive, B° otherwise. 

If both lepton and kaon tags are present and 
provide inconsistent flavor tags, the event is re- 
jected from the lepton and kaon tagging cate- 
gories. 
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5.2 Neural network categories 

The use of a second tagging algorithm is moti- 
vated by the potential flavor-tagging power car- 
ried by non-identified leptons and kaons, correla- 
tions between leptons and kaons, multiple kaons, 
softer leptons from charm semileptonic decays, 
soft pions from D* decays and more generally 
by the momentum spectrum of charged parti- 
cles from B meson decays. One way to exploit 
the information contained in a set of correlated 
quantities is to use multivariate methods such 
as neural networks. 

We define five different neural networks, 
called feature nets, each with a specific goal. 
Four of the five feature nets are track-based : 
the L and LS feature nets are sensitive to the 
presence of primary and cascade leptons, respec- 
tively, the K feature net to that of charged kaons 
and the Sof tPi feature net to that of soft pions 
from D* decays. In addition, the Q feature net 
exploits the charge of the fastest particles in the 
event. 

The variables used as input to the neural net- 
work tagger are the highest values of the L, LS 
and SoftPi feature net outputs multiplied by 
the charge, the highest and the second highest 
value of the K feature net output multiplied by 
the charge, and the output of the Q feature net. 

The output of the neural network tagger, 
xnt, can be mapped onto the interval [—1,1]. 
The tag is B° if xnt is negative, B° other- 
wise. Events with |xtvt| > 0.5 are classified 
in the NT1 tagging category and events with 
0.2 < \xnt\ < 0.5 in the NT2 tagging category. 
Events with I^ntI < 0.2 have very little tagging 
power and are excluded from the sample used in 
the analysis. 

6 Measurement of mistag frac- 
tions 



The mistag fractions are measured directly in 
events in which one B° candidate, called the 
B rec , is fully reconstructed in a flavor eigenstate 
mode. The flavor-tagging algorithms described 



in the previous section are applied to the rest of 
the event, which constitutes the potential Bt ag . 

Considering the B°B° system whole, one 
can classify the tagged events as mixed or un- 
mixed depending on whether the B tag is tagged 
with the same flavor as the B rec or with the 
opposite flavor. Neglecting the effect of pos- 
sible background contributions, and assuming 
the B rec is properly tagged, one can express 
the measured time-integrated fraction of mixed 
events \ as a function of the precisely-measured 
B°B° mixing probability Xd ■ 



X = Xd + (1 - 2xd) w 



(8) 



where Xd = \ x 2 d /(l + x 2 d ), with x d = Am d /T. 
Thus one can deduce an experimental value of 
the mistag fraction w from the data. 

A time-dependent analysis of the fraction of 
mixed events is even more sensitive to the mistag 
fraction. The mixing probability is smallest for 
small values of At = t rec — tf ag so that the ap- 
parent rate of mixed events near At=0 is gov- 
erned by the mistag probability (see Fig. [6j). A 
time-dependent analysis can also help discrimi- 
nate against backgrounds with different time de- 
pendence. 

By analogy with Eq. |2|, we can express the 
density functions for unmixed (+) and mixed 
(— ) events as 



H±{At; T, Am d , V, a) = 

h±(At; T, Am d , V)®K( At; 



o) 



0) 



where 



h±(At; T, Am d , V) = 

re _r|At| [1 ± P x cosAm d Ai]. (10) 

These functions are used to build the log- 
likelihood function for the mixing analysis: 



ln£ 



M 



Eunmixed ln ^+ ( *> ^ \ Am d , a, T>i 
+ Emixed ln ^-U; r > Am d, a, Vi) 
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Figure 6: Fraction of mixed events m/(u + m) 
as a function of |At| (ps) for data events in 
the hadronic sample, for neutral B mesons (full 
squares) and charged B mesons (open circles). 
All tagging categories are included. This rate 
is a constant as a function of At for charged 
B mesons, but develops a mixing oscillation for 
neutral B mesons. The rate of mixed events 
extrapolated to At = is governed by the 
mistag fraction w. The dot-dashed line at t cu t = 
2.5 ps indicates the bin boundary of the time- 
integrated single-bin method. 



which is maximized to extract the estimates of 
the mistag fractions Wi = ^(1 — T>i). 

The extraction of the mistag probabilities for 
each tagging category is complicated by the pos- 
sible presence of mode-dependent backgrounds. 
We deal with these by adding specific terms in 
the likelihood functions describing the different 
types of backgrounds (zero lifetime, non-zero 
lifetime without mixing, non-zero lifetime with 
mixing), described in Ref. 6. 

A simple time-integrated single bin method 
is used ELS cL check of the time-dependent anal- 
ysis for the determination of dilutions from the 
fully reconstructed B° sample. The mistag frac- 
tions are deduced from the number of unmixed 
events, u, and the number of mixed events, m, 
in a single optimized At interval, | At | < t cu t- 



The bin boundary t cu t, chosen to minimize the 
statistical uncertainty on the measurement, is 
equal to 2.5 ps, i.e., 1.6 B° lifetimes. (t C ut is 
indicated by a dot-dashed line in Fig. |6[) The 
resulting mistag fractions based on this method 
are in good agreement with the mistag fractions 
obtained with the maximum-likelihood fit||. 

6.1 Tagging efficiencies and mistag 
fractions 

The mistag fractions and the tagging efficiencies 
are summarized in Table ||. We find a tagging 
efficiency of (76.7±0.5)% (statistical error only). 
The lepton categories have the lowest mistag 
fractions, but also have low efficiency. The Kaon 
category, despite having a larger mistag fraction 
(19.7%), has a higher effective tagging efficiency; 
one-third of events are assigned to this category. 
Altogether, lepton and kaon categories have an 
effective tagging efficiency Q ~ 20.8%. Most of 
the separation into B° and B° in the NT1 and 
NT2 tagging categories derives from the Sof tPi 
and Q feature nets. Simulation studies indi- 
cate that roughly 40% of the effective tagging 
efficiency occurs in events that contain a soft 
7r aligned with the Bt ag thrust axis, 25% from 
events which have a track with p* > 1.1 GeV/c, 
10% from events which contain multiple leptons 
or kaons with opposite charges and are thus not 
previously used in tagging, and the remaining 
25% from a mixture of the various feature nets. 
The neural network categories increase the ef- 
fective tagging efficiency by ~ 7% to an overall 
Q = (27.9 ± 1.6)% (statistical error only). 

Of the 168 CP candidates, 120 are tagged: 
70 as B° and 50 as B°. The number of tagged 
events per category is given in Table [|. 

7 Systematic uncertainties and 
cross checks 

Systematic errors arise from uncertainties in 
input parameters to the maximum likelihood 
fit, incomplete knowledge of the time resolu- 
tion function, uncertainties in the mistag frac- 
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Table 3: Mistag fractions measured from a maximum-likelihood fit to the time distribution for the 
fully-reconstructed B° sample. The Electron and Muon categories are grouped into one Lepton 
category. The uncertainties on e and Q are statistical only. 



Tagging Category 


£(%) 


w (%) 


Q{%) 


Lepton 
Kaon 
NT1 
NT2 


11.2 ±0.5 
36.7 ±0.9 
11.7 ±0.5 
16.6 ±0.6 


9.6 ± 1.7 ± 1.3 
19.7 ± 1.3 ± 1.1 
16.7 ±2.2 ±2.0 
33.1 ±2.1 ±2.1 


7.3 ±0.7 
13.5 ± 1.2 
5.2 ±0.7 
1.9 ±0.5 


all 


76.7 ±0.5 




27.9 ± 1.6 



Table 4: Categories of tagged events in the CP sample. 









J/i> 








t/j(2S)K° s 


CP sample 


Tagging Category 


(K° s 


— ► 7T~ 


f 7T-) 


(K°s 


- ttV) 




— ► 7T + 7T ) 


(tagge 


d) 




B° 


B° 


all 


B° 


B° 


all 


B° 


B° 


all 


B° 


B° 


all 


Electron 


1 


3 


4 


1 





1 


1 


2 


3 


3 


5 


8 


Muon 


1 


3 


4 











2 





2 


3 


3 


6 


Kaon 


29 


18 


47 


2 


2 


4 


5 


7 


12 


36 


27 


63 


NT1 


9 


2 


11 


1 





1 


2 





2 


12 


2 


14 


NT2 


10 


9 


19 


3 


3 


6 


3 


1 


4 


16 


13 


29 


Total 


50 


35 


85 


7 


5 


12 


13 


10 


23 


70 


50 


120 



tions, and possible limitations in the analysis 
procedure. We fix the B° lifetime to the nom- 
inal PDGjlCfl central value Tgo = 1.548 ps and 
the value of Am^ to the nominal PDG value 
Arrid = 0.472 ?ips _1 . The errors on sin2/3 due 
to uncertainties in Tgo and Arrid are 0.002 and 
0.015, respectively. The remaining systematic 
uncertainties are discussed in the following sec- 
tions. 

7.1 Uncertainties in the resolution 
function 

The time resolution is measured with the 
high-statistics sample of fully-reconstructed B° 
events. The time resolution for the CP sample 
should be very similar, especially to that mea- 
sured for the hadronic sample. We verify that 
the resolution function extracted in the hadronic 



sample is consistent with the one extracted in 
the semileptonic sample. We assign as a sys- 
tematic error the variation in sin2/3 obtained by 
changing the resolution parameters by one sta- 
tistical standard deviation. The corresponding 
error on sin2/3 is 0.019. 

We use a full Monte Carlo simulation to ver- 
ify that the Bremsstrahlung recovery procedure 
in the J/tp — > e + e~ mode does not introduce 
any systematic bias in the At measurement, nor 
does it affect the vertex resolution and pull dis- 
tributions. 

In order to check the impact of imperfect 
knowledge of the bias in At on the measure- 
ment, we allow the bias of the second Gaussian 
to increase to 0.5 ps. The resulting change in 
sin2/3 of 0.047 is assigned as a systematic error. 
The sensitivity to the bias is due to the different 



12 



number of events tagged as B° and B° . 

7.2 Uncertainties in flavor tagging 

The mistag fractions are measured with uncer- 
tainties that are either correlated or uncor re- 
lated between tagging categories. We study the 
effect of uncorrelated errors (including statis- 
tical errors) on the asymmetry by varying the 
mistag fractions individually for each category, 
using the full covariance matrix. For correlated 
errors, we vary the mistag fractions for all cate- 
gories simultaneously. 

The main common source of systematic un- 
certainties in the measurement of mistag frac- 
tions is the presence of backgrounds, which are 
more significant in the semileptonic sample than 
in the hadronic sample. The largest background 
is due to random combinations of particles and 
can be studied with mass sidebands. Additional 
backgrounds arise in the semileptonic sample 
from misidentified leptons, from leptons incor- 
rectly associated with a true D* from B decays, 
and from charm events containing a D* and a 
lepton. The details of the procedure for account- 
ing for the backgrounds and the uncertainties on 
the background levels, and the estimates of re- 
sulting systematic errors on the mistag fractions 
are given in Ref. 6. We estimate the systematic 
error on sin2/3 due to the uncertainties in the 
measurement of the mistag fractions to be 0.053, 
for our CP sample. 

In the likelihood function, we use the same 
mistag fractions for the B° and B° samples. 
However, differences are expected due to effects 
such as the different cross sections for K + and 
K~ hadronic interactions. For equal numbers of 
tagged B° and B° events, the impact on sin2/3 of 
a difference in mistag fraction, 5w = w B o—w^ , is 
insignificant. From studies of charged and neu- 
tral B samples, we find that the mistag differ- 
ences are < 0.02 for the NT1 category, < 0.04 for 
the Kaon category, and negligible for the lepton 
categories. However, for the NT2 category, there 
is a significant difference between the B° and B° 
mistag fractions, 5w = 0.16, which is not pre- 
dicted by our simulation. Although this would 



lead to a negligible systematic shift in sin2/?, we 
cover the possibility of different mistag fractions 
in the CP sample and the fully-reconstructed 
sample used to measure the mistag fractions by 
assigning as a systematic uncertainty the shift in 
sin2/3 resulting from using the measured mistag 
fraction for the NT2 category from the sample of 
J/i/jK*° events only. The resulting conservative 
systematic uncertainty on sin2/3 is 0.050. 

For a small sample of events, there can be a 
significant difference in the number of B° and 5° 
events, AN = N B o — N^ . For a single tagging 
category, the fractional change in sin2/3 from 
such a difference is Asin2/3/ sin2/3 ~ 5wAN/N . 
In the CP sample, AN/N is significant only in 
the Kaon and NT1 categories (see Table |]). Tak- 
ing into account their relative weight in the over- 
all result, we assign a fractional systematic error 
of 0.005 on sin2/3. 

The systematic uncertainties on the mistag 
fractions due to the uncertainties on t b o and 
Am,i are negligible. 

7.3 Uncertainties due to backgrounds 

The fraction of background events in the CP 
sample (J/ip K® and ip(2S)Kg) is estimated to 
be (5±3)%. The portion of this background that 
occurs at small values of At (e.g., contributions 
from u, d and s continuum events) does not con- 
tribute substantially to the determination of the 
asymmetry. We estimate that this reduces the 
effective background to 3%. We correct for the 
background by increasing the apparent asymme- 
try by a factor of 1.03. In addition, we assign a 
fractional systematic uncertainty of 3% on the 
asymmetry, to cover both the uncertainty in the 
size of the background and the possibility that 
the background might have some CP-violating 
component. 

8 Extracting sin2/3 

8.0.1 Blind analysis 

We have adopted a blind analysis for the ex- 
traction of sin2/3 in order to eliminate possible 
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experimenter's bias. We use a technique that 
hides not only the result of the unbinned maxi- 
mum likelihood fit, but also the visual CP asym- 
metry in the At distribution. The error on the 
asymmetry is not hidden. 
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Figure 7: Variation of the log likelihood 
function of sin2/3. The two horizontal dashed 
lines indicate changes in the log likelihood cor- 
responding to one and two statistical standard 
deviations. 

The amplitude of the asymmetry Acp(At) 
from the fit is hidden from the experimenter by 
arbitrarily flipping its sign and adding an ar- 
bitrary offset. The sign flip hides whether a 
change in the analysis increases or decreases the 
resulting asymmetry. However, the magnitude 
of the change is not hidden. 

The visual CP asymmetry in the At distri- 
bution is hidden by multiplying At by the sign 
of the tag and adding an arbitrary offset. 

With these techniques, systematic studies 
can be performed while keeping the numerical 
value of sin2/3 hidden. In particular, we can 
check that the hidden At distributions are con- 
sistent for B and B° tagged events. The same is 
true for all the other checks concerning tagging, 
vertex resolution and the correlations between 
them. For instance, fit results in the differ- 
ent tagging categories can be compared to each 
other, since each fit is hidden in the same way. 
The analysis procedure for extracting sin2/3 was 
frozen, and the data sample fixed, prior to un- 
blinding. 
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Figure 8: Distribution of At for (a) the B° 
tagged events and (b) the B° tagged events 
in the CP sample. The error bars plotted for 
each data point assume Poisson statistics. The 
curves correspond to the result of the unbinned 
maximum-likelihood fit and are each normal- 
ized to the observed number of tagged B° or 
B° events. 



8.1 Cross checks of the fitting proce- 
dure 

We submitted our maximum-likelihood fitting 
procedure to an extensive series of simulation 
tests. The tests were carried out with two differ- 
ent implementations of the fitting algorithm to 
check for software errors. The validation stud- 
ies were done on two types of simulated event 
samples. 

• "Toy" Monte Carlo simulation tests. In 
these samples, the detector response is not 
simulated. Monte Carlo techniques are 
used with parametrized resolution func- 
tions and tagging probabilities. We vali- 
dated the fitting procedure on large sam- 
ples of simulated CP events, for various 
numbers of tagging categories, values of 
mistag fractions and values of sin2/3. We 
also simulated a large number of 100-event 
experiments, with the purpose of investi- 
gating statistical issues with small sam- 
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Table 5: Result of fitting for CP asymmetries in 
the entire CP sample and in various subsamples. 



At (ps) 

Figure 9: The raw B°-B° asymmetry (Ngo — 
Ng )/(N B o + iVgo), with binomial errors, is 
shown as a function of At. The time-dependent 
asymmetry is represented by a solid curve for 
our central value of sin2/3, and by two dotted 
curves for the values at plus and minus one 
statistical standard deviation from the central 
value. The curves are not centered at (0, 0) in 
part because the probability density functions 
are normalized separately for B° and B° events, 
and our CP sample contains an unequal num- 
ber of B° and B° tagged events (70 B° versus 
50 B°). The x 2 between the binned asymmetry 
and the result of the maximum-likelihood fit is 
9.2 for 7 degrees of freedom. 



pies, including values of sin2/3 near un- 
physical regions. We checked that the fit- 
ter performs well in the presence of back- 
grounds for the extraction of the mistag 
fractions. We exercised the combined CP 
and mixing fits, and found that although 
combined fits perform well, they do not 
significantly improve the statistical sensi- 
tivity of the result. 

Full Monte Carlo simulation tests. We 
studied samples of J/ipK®, J/tpK + , D*ir 
and D*£v events produced with the BABAR 
GEANT3 detector simulation and recon- 



sample 


smzp 


CP sample 


0.12±0.37 


J/^K% (K% -> vr+vr^) 


-0.10 ±0.42 


other CP events 


0.87 ±0.81 


Lepton 


1.6 ± 1.0 


Kaon 


0.14 ±0.47 


NT1 


-0.59 ±0.87 


NT2 


-0.96 ± 1.30 



structed with the BABAR reconstruction 
program. J/tpK®, events were generated 
with various values of sin2/3. We ex- 
tracted the "apparent CP-asymmetry" for 
the charged -B's and found it to be con- 
sistent with zero. We studied the differ- 
ence in tagging efficiencies and in mistag 
fractions between the charged and neutral 
B samples. We also tested the procedure 
for extracting the mistag fractions from 
hadronic and semileptonic samples of fully 
simulated events (D*tt and D*£v). 

9 Results 

The maximum-likelihood fit for sin2/3, using the 
full tagged sample of B°/B° — ► J/ipK® and 
B°/B° -» 4>{2S)K® events, gives: 

sin2/3 = 0.12 ± 0.37(stat) ± 0.09(syst) . 

For this result, the B° lifetime and Am<j are 



fixed to the current best values [10], and At reso- 
lution parameters and the mistag rates are fixed 
to the values obtained from data as summarized 
in Tables |2| and ||. The log likelihood is shown 
as a function of sin2/3 in Fig. ^, the At distribu- 
tions for B° and B° tags in Fig. || and the raw 
asymmetry as a function of At in Fig. |9[ The 
results of the fit for each type of CP sample and 
for each tagging category are given in Table |5[ 
The contributions to the systematic uncertainty 
are summarized in Table 0. 
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Table 6: Summary of systematic uncertainties. We compute the fractional systematic errors 
using the actual value of our asymmetry increased by one statistical standard deviation, that is 
0.12 + 0.37 = 0.49. The different contributions to the systematic error are added in quadrature. 



source oi uncertainty 


Uncertainty on sin2/? 


TT i ■ j I) 

Uncertainty on Tg 
Uncertainty on Am^ 

Uncertainty on Az resolution for CP sample 
Uncertainty on time-resolution bias for CP sample 
Uncertainty on measurement of mistag fractions 
Different mistag fractions for CP and non-CP samples 
Different mistag fractions for B° and B° 
Background in CP sample 


0.002 
0.015 
0.019 
0.047 
0.053 
0.050 
0.005 
0.015 


Total systematic error 


0.091 



We estimate the probability of obtaining 
the observed value of the statistical uncertainty, 
0.37, on our measurement of sin2/3 by generating 
a large number of toy Monte Carlo experiments 
with the same number of tagged CP events, and 
distributed in the same tagging categories, as in 
the CP sample in the data. We find that the 
errors are distributed around 0.32 with a stan- 
dard deviation of 0.03, and that the probability 
of obtaining a value of the statistical error larger 
than the one we observe is 5%. Based on a large 
number of full Monte Carlo simulated experi- 
ments with the same number of events as our 
data sample, we estimate that the probability 
of finding a lower value of the likelihood than 
our observed value is 20%. 

10 Validating analyses 

To validate the analysis we use the charmonium 
control sample, composed of B + — > J/ipK + 
events and events with self-tagged J/ip K*° 
(K*° -> K+tt-) neutral P's. We also use 
the event samples with fully-reconstructed can- 
didates in charged or neutral hadronic modes. 
These samples should exhibit no time-dependent 
asymmetry. In order to investigate this exper- 
imentally, we define an "apparent CP asymme- 
try" , analogous to sin2/3 in Eq. ||, which we ex- 



Table 7: Results of fitting for apparent 
CP asymmetries in various charged or neutral 
flavor-eigenstate B samples. 



Sample 


Apparent 
CP-asymmetry 


Hadronic B ± decays 


0.03 ± 0.07 


Hadronic B° decays 


-0.01 ±0.08 


J/ipK+ 


0.13 ±0.14 


K*° _> K + TT~ 


0.49 ± 0.26 



tract from the data using an identical maximum- 
likelihood procedure. 

The events in the control samples are flavor 
eigenstates and not CP eigenstates. They are 
used for testing the fitting procedure with the 
same tagging algorithm as for the CP sample 
and, in the case of the B + modes, with self- 
tagging based on their charge. We also perform 
the fits for B° and B° (or B + and B~) events 
separately to study possible flavor-dependent 
systematic effects. For the charged B modes, 
we use mistag fractions measured from the sam- 
ple of hadronic charged B decays. 

In all fits, including the fits to charged sam- 
ples, we fix the lifetime r^o and the oscillation 
frequency Amj to the PDG values [|lC|] . 
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The results of a series of validation checks on 
the control samples are summarized in Table fj]. 

The two high-statistics samples and the 
J/ipK + sample give an apparent CP asymme- 
try consistent with zero. The 1.9 a asymmetry 
in the J/ip K*° is interpreted as a statistical fluc- 
tuation. 

Other BABAR time-dependent analyses pre- 
sented at this Conference demonstrate the valid- 
ity of the novel technique developed for use at an 
asymmetric B Factory. The measurement of the 
B°-B° oscillation frequency described in Ref. 6 
uses the same time resolution function and tag- 
ging algorithm as the CP analysis. Fitting for 
Anid in the maximum-likelihood fit for the fully- 
reconstructed hadronic and semileptonic neutral 
B decays, we measure 

Am d = 

0.512 ± 0.017(stat) ± 0.022(syst) h ps" 1 , 

which is consistent with the world average [[lOf| 
Am d = 0.472 ± 0.017 Tips" 1 . The B° lifetime 
measurement described in Ref. 7 uses the same 
inclusive vertex reconstruction technique as the 
CP analysis. We measure 

t b o = 1.506 ± 0.052(stat) ± 0.029(syst) ps , 

also consistent with the world average t b o = 
1.548 ± 0.032 ps. 

11 Conclusions and prospects 

We have presented BABAR' 1 s first measurement 
of the CP-violating asymmetry parameter sin2/3 
in the B meson system: 

sin2/3 = 0.12 ± 0.37(stat) ± 0.09(syst) . 

Our measurement is consistent with the world 
average^] sin2/3 = 0.9 ± 0.4, and is currently lim- 
ited by the size of the CP sample. We expect 
to more than double the present data sample in 
the near future. 

3 Based on the OPAL result jn) sin2/3 = 3.2±|^ ± 0.5 
and the CDF resultjt| sin2fl = Q.TS^m- See also 
ALEPH's preliminary result @ sin2/3 = 0.93±^ || ± 1|. 
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Figure 10: Present constraints on the posi- 
tion of the apex of the Unitarity Triangle in the 
(p, fj) plane. The fitting procedure is described 
in Ref. 2. Our result sin2/3 = 0.12 ± 0.37(stat) 
is represented by cross-hatched regions corre- 
sponding to one and two statistical standard de- 
viations. 



Figure [T(] shows the Unitarity Triangle in 
the (p, fj) plane, with BABAR 's measured central 
value of sin2/3 shown as two straight lines; there 
is a two-fold ambiguity in deriving a value of (3 
from a measurement of sin2/3 . Both choices are 
shown with cross-hatched regions corresponding 
to one and two times the one-standard-deviation 
experimental uncertainty. The ellipses corre- 
spond to the regions allowed by all other mea- 
surements that constrain the Unitarity Trian- 
gle. Rather than make the common, albeit un- 
founded, assumption that our lack of knowledge 
of theoretical quantities, or differences between 
theoretical models, can be parametrized (typi- 
cally Gaussian or flat distribution), we have 
chosen to display the ellipses corresponding to 
measurement errors at a variety of representa- 
tive choices^ of theoretical parameters. The fit- 

4 We use the following set of measurements: \V c b\ = 
0.0402 ± 0.017, \V ub /V cb \ = (\Vub/Vcb\) ± 0.0079, Am d = 
0.472 ± 0.017 ftps -1 and \e K \ = (2.271 ± 0.017) x 
10 -3 , and for Am s the set of amplitudes correspond- 
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ting procedure is described in Ref. 2. 

While the current experimental uncertainty 
on sin2/3 is large, the next few years will bring 
substantial improvements in precision, as well as 
measurements for other final states in which CP- 
violating asymmetries are proportional to sin2/3, 
and measurements for modes in which the asym- 
metry is proportional to sin2<x 
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ing to a 95%CL limit of 14.6 ftps" 1 . We scan the 
model-dependent parameters (| V u b/V c b\), Bk, fB d \/BB d 
and f s , in the range [0.070, 0.100 ], [0.720, 0.980 ], 
[185, 255] MeV and [1.07, 1.21], respectively. 



References 

[1] J. H. Christenson et al, Phys. Rev. Lett. 
13, 138 (1964); 

NA31 Collaboration, G. D. Barr et al, 
Phys. Lett. 317, 233 (1993); 
E731 Collaboration, L. K. Gibbons et al, 
Phys. Rev. Lett. 70, 1203 (1993). 

[2] See, for instance, "Overall determinations 
of the CKM matrix", in P. H. Harrison 
and H. R. Quinn, eds., "The BABAR physics 
book", SLAC-R-504 (1998) section 14 and 
references therein. 

[3] For an introduction to CP violation, see, 
for instance, "A CP violation primer", in 
P. H. Harrison and H. R. Quinn, eds., "The 
BABAR physics book", SLAC-R-504 (1998) 
section 1 and references therein. 

[4] C. Jarlskog, in CP Violation, C. Jarlskog 
ed., World Scientific, Singapore (1988). 

[5] BABAR Collaboration, B. Aubert et al, 
"Exclusive B decays to charmonium fi- 
nal states", BABAR-CONF-00/05, submit- 
ted to the XXX th International Conference 
on High Energy Physics, Osaka, Japan. 

[6] BABAR Collaboration, B. Aubert et al, "A 
measurement of the B°B° oscillation fre- 
quency and determination of flavor-tagging 
efficiency using semileptonic and hadronic 
B decays", BABAR-CONF-00/08, submit- 
ted to the XXX th International Conference 
on High Energy Physics, Osaka, Japan. 

[7] BABAR Collaboration, B. Aubert et al, "A 
measurement of the charged and neutral 
B meson lifetimes using fully reconstructed 
decays", BABAR-CONF-00/07, submitted 
to the XXX th International Conference on 
High Energy Physics, Osaka, Japan. 

[8] BABAR Collaboration, B. Aubert et al, 
"The first year of the BABAR experiment 
at PEP-II", BABAR-CONF-00/17, submit- 
ted to the XXX th International Conference 
on High Energy Physics, Osaka, Japan. 



18 



[9] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 
41, 1581 (1978). 

[10] Particle Data Group, D. E. Groom et ai, 
Eur. Phys. Jour. C 15, 1 (2000). 

[11] OPAL Collaboration, K. Ackerstaff et ai, 
Eur. Phys. Jour. C 5, 379 (1998). 

[12] CDF Collaboration, T. Affolder et al, 
Phys. Rev. D61, 072005 (2000). 

[13] ALEPH Collaboration, ALEPH 99-099 
CONF 99-54 (1999). 



19 



